1. Introduction 

The search for the top quark and the Higgs boson, the only as yet unobserved ingre- 
dients of the Standard Model (SM), is expected to dominate the area of particle physics 
in the coming years. While the Higgs boson has always been naturally associated to the 
breaking of the electro-weak symmetry, the recent experimental discovery that the top 
quark is heavier than the W and Z has led physicists to speculate that it also may be sen- 
sitive to physics on the electro-weak scale. For this reason, the theoretical determination 
of its physical properties in the context of the SM and their systematic comparison with 
experiment may reveal an underlying, more fundamental theory. 

One of the most efficient ways to study top quarks will be to pair-produce them in 
future very energetic e + e~ colliders [1], through the reaction e + e~ — > ti. In general, 
the leptonic nature of the target allows for clean signals. In addition, due to their large 
masses, the produced top quarks are expected to decay weakly (ti — > bW + bW~ , with 
subsequent leptonic decays of the W), before hadronization takes place; therefore electro- 
weak properties of the top can be studied in detail and QCD corrections can be reliably 
evaluated in the context of perturbation theory, when the energy of the collider is well 
above the threshold for ti production. 

Motivated by this type of experiments, Atwood and Soni presented in a recent paper 
[2] a phenomenological analysis for determining the magnetic and electric dipole moment 
form factors of the top quark. Such form factors are defined through the following Lorentz 
decomposition of the Vti vertex, where V represents a boson (a 7 or Z in our case) coupled 
to the conserved leptonic current: 

r^ 2 ) = 7„W) + *SV) + 7m76*JV) + ^ 7**2V) (1-1) 

where a^ v = | [7^,7^] and q 2 = s, i.e. the Mandelstam variable associated to the squared 
energy of the center of mass. In the above decomposition, F% is the magnetic dipole 
moment (MDM) and FY is the electric dipole moment (EDM) form factor. In particular, 
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F^ defined at q 2 = is the usual definition of the anomalous magnetic moment. In the case 
of the top quark production, clearly q 2 > Am 2 . Within the SM the tree- level value for both 
F2 and FY is zero. F^ becomes non-zero through one-loop quantum corrections, whereas 
F± , which violates CP, receives its first non-vanishing contributions at three loops [3]. 
The upshot of the analysis of [2] was that the dependence of the differential cross section 
for the reaction e + e~ — > tt on the real and imaginary (absorptive) parts of the MDM and 
EDM form factors, for an incoming photon or Z, can be singled out individually, through 
a set of optimally chosen physical observables. Clearly, the possibility of such a detailed 
experimental study of the top dynamics is attractive, and could provide valuable probes 
for new physics. 

The theoretical aspects of the situation are however not entirely clear. As it was 
pointed out already in the classic paper by Fujikawa, Lee, and Sanda [4], off-shell form 
factors of fermions are in general gauge dependent quantities. In the context of the i?g 
gauges, for example, a residual dependence on the gauge-fixing parameter £ survives in 
the final expressions of form factors, when q 2 7^ [5]. Obviously, in the case of e + e~ 
annihilation into heavy fermions, the value of q 2 must be above the heavy fermion threshold 
(q 2 > Am 2 , in our case). Consequently, the intermediate photon and Z are far off-shell, 
and therefore, MDM and EDM form factors may in general be gauge-dependent and not 
suitable for comparison with experiment. 

A way out of this problem is to define gauge independent off-shell vertices using the 
pinch technique (PT). The PT was invented by Cornwall over a decade ago [6] and has 
since been applied to a variety of physical problems. The main idea of this method is to 
resum via a well-defined algorithm the Feynman diagrams contributing to an ostensibly 
gauge-invariant process (like an S- matrix element), in such a way as to form new gauge- 
independent proper vertices, and new propagators with gauge-independent self-energies 
and only a trivial gauge dependence - that of their tree level parts. In the context of 
QCD a gauge invariant gluon self-energy was derived, and its Schwinger-Dyson equation 
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constructed and solved for T = [7] , as well as finite T [8] . The plasmon decay rate was 
also calculated at finite T using the same method [9]. Later the QCD gauge invariant 
three-gluon vertex was calculated at one loop level and was shown to satisfy a simple QED- 
like Ward identity [10]. The subleading corrections to the self-energy were calculated by 
Lavelle [11]. Finally, the gauge-invariant four-gluon QCD vertex was constructed and 
its Ward identity derived [12]. The PT was first extended to the case of non-Abelian 
gauge theories with spontaneously broken gauge symmetry (with elementary Higgs) in the 
context of a toy field theory based on SU(2), and a gauge independent electro-magnetic 
form factor for the neutrino was constructed [13]. The complicated task of applying 
the P.T. in the electro-weak sector of the Standard Model was accomplished by Degrassi 
and Sirlin [14]. These last authors, in addition to deriving explicit expressions for the 
one loop gauge-invariant WW and ZZ self-energies, introduced an alternative description 
of the PT in terms of equal time commutators of currents. More recently [15], gauge 
independent *yW + W~ and ZW + W~ vertices were constructed in the context of the SM; 
they are related to the gauge-independent W-propagators introduced in [14] by a very 
simple Ward identity, and give rise to expressions for the W magnetic dipole and electric 
quadrupole moments, which, unlike previous treatments [16], are independent of the 
gauge-fixing parameter and infrared finite. 

In this paper we focus on the MDM form factor of the top quark. In particular: 
a. In the context of the SM we show that the conventionally defined form factor for the 
top quark is gauge dependent in the class of gauges. The gauge dependence is formally 
expressed as a double integral over the usual Feynman parameters, and is proportional to 



b. We show how one can define a gauge independent MDM form factor using the PT. 

c. The gauge dependence found in (a) is computed and turns out to be numerically 
very strong; its presence distorts not only the quantitative but also the qualitative behavior 
of the answer obtained in (b). More specifically, unphysical thresholds are introduced, 

3 



the asymptotic behavior of the form factors as q 2 — > oo is altered, and the numerical 
dominance of perturbative QCD, which is present in the gauge independent expressions 
of (b), is totally washed out. Of particular interest is the fact that the popular "unitary" 
gauge ( the limit of the gauges as £ — > oo) gives a completely wrong answer. This 
analysis indicates in retrospect that the gauge dependence established in (a) is a serious 
pathology and may lead to erroneous conclusions. 

d. We extract the real and imaginary parts of the gauge independent MDM form 
factors for incoming photon and Z (i 7 ^ and , respectively), and evaluate them numeri- 
cally. We perform our computations for values of q 2 comfortably above the threshold for 
top production, and we ignore therefore non-perturbative effects due to threshold dynamics 
[17]. It turns out that even though perturbative QCD provides in general the dominant 
contribution, the MDM form factor displays a high sensitivity to the Higgs mass. It is 
important to emphasize that the Higgs dependence of the MDM is actually stronger than 
that displayed by the e + e~ — > ti cross section as a whole [18]. 

The paper is organized as follows. In section 2 we review the S-matrix PT and discuss 
some of the more important results for our purposes. In addition, we briefly present 
Degrassi's and Sirlin's alternative formulation of the PT. In section 3 we perform an explicit 
calculation and determine the naive, gauge dependent expression for the MDM form factor. 
In section 4 we show how to use the pinch technique to define gauge independent r yti and 
Zti vertices, which give rise to gauge independent expressions for the MDM form factors. 
In section 5 we present a simple QED Ward identity relating the gauge independent jti 
vertex with the inverse top quark propagators, calculated in the Feynman gauge (£ = 1). 
Section 6 contains the results of our gauge invariant MDM computations, and a discussion 
of their significance. We present our conclusions in section 7. Finally, in an Appendix, we 
list in detail the contributions of individual Feynman diagrams to the MDM form factors. 
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2. The Pinch Technique 

In this section we briefly review the S-matrix pinch technique. In particular, we 
outline the method of derivation of the gauge-independent proper self-energy of a gauge 
boson and comment on the technical differences that arise when PT is applied to a theory 
with symmetry breaking, like the electro-weak theory, as opposed to a theory like QCD. 
In addition, we present the main idea of Degrassi's and Sirlin's formulation of the pinch, 
and establish some of the notation we will use in the sequel. 

The PT is an algorithm that allows the construction of modified gauge independent n- 
point functions, through the order by order resummation of Feynman graphs contributing 
to a certain physical and therefore ostensibly gauge independent process (an S-matrix 
in our case). The simplest example that demonstrates how the P.T. works is the gauge 
boson two point function (propagator). Consider the S'-matrix element T for the elastic 
scattering of two fermions of masses Mi and M 2 . To any order in perturbation theory T 
is independent of the gauge fixing parameter £, defined by the free gluon propagator 

Aju/Cg) = - 2 — (2.1) 

On the other hand, as an explicit calculation shows, the conventionally defined proper 
self-energy (collectively depicted in graph la) depends on £. At the one loop level this 
dependence is cancelled by contributions from other graphs, like lb and lc, which do not 
seem to be of propagator type at first glance. That this must be so is evident from the 
form of T: 

T(s, t, Mi, M 2 ) = Tx (£) + T 2 (t, M u M 2 ) + T 3 (s, t, M u M 2 ) (2.2) 

where the function Ti(t) depends only on the Mandelstam variable t = — (pi—pi) 2 = —q 2 , 
and not on s = (pi + p 2 ) 2 or on the external masses. The propagator-like parts of graphs 
like le and If, which enforce the gauge independence of Ti(t), are called "pinch parts". 
The pinch parts emerge every time a gluon propagator or an elementary three-gluon vertex 
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contribute a longitudinal k^ to the graph's numerator. The action of such a term is to 
trigger an elementary Ward identity of the form 



once it gets contracted with a 7 matrix. The first term on the right-hand side of Eq. (2.3) 
will remove the internal fermion propagator - that is a "pinch" - whereas S~ l (p) vanishes 
on shell. This last property characterizes the S-matrix PT we will use throughout this 
paper. Returning to the decomposition of Eq. (2.2), the function T\(t) is gauge invariant 
and unique, and represents the contribution of the new propagator. We can construct the 
new propagator, or equivalently Ti(t), directly from the Feynman rules. In doing so it is 
evident that any value for the gauge parameter £ may be chosen, since Ti, T 2 , and T 3 are 
all independent of £. The simplest of all covariant gauges is certainly the Feynman gauge 
(£ = 1), which gets rid of the longitudinal part of the gluon propagator. Therefore, the only 
possibility for pinching in this gauge arises from the four-momentum of the three-gluon 
vertices, and the only propagator-like contributions come from graph lb. 

To explicitly calculate the pinching contribution of a graph such as lb it is convenient 
to decompose the vertex in the following way, first proposed by 't Hooft. Group theory 
factors aside, 



k"-y li =t = y + t-m)-y-m) 



(2.3) 



S-\p + k) -s-\p) 



r — r p -4- r F 



(2.4) 



with 




(q + k) v g lia + k^g, 



(2.5) 



and 




(2.6) 



r p a satisfies a Feynman-gauge Ward identity: 




[k 2 ~ (k + qf]g^ 



(2.7) 
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where the RHS is the difference of two inverse propagators in the Feynman gauge. As for 
T^ va (P for "pinch"), it gives rise to pinch parts when contracted with 7 matrices 

9md + V) =9n*[(j> + 4-m)-ty-$-m)] 

(2.8) 

=9va[S- 1 (p + q)-S- 1 (p-k)] 

and 

9va# =9va[(]f> - m) - - $ - m)] 

(2.9) 

=9ua[S- 1 ( P )-S- 1 (p-k)] 

Now both S~ 1 (p + q) and S~ 1 (p) vanish on shell, whereas the two terms proportional to 
S~ 1 (p — k) pinch out the internal fermion propagator in graph lb, and so we are left with 
two "pinch" (propagator-like) parts and one "regular" (purely vertex-like) part: 

pinch part = 2 x (_I C>1 )[^] J fc2( f* g)2 (<7 W 7„) (2.10) 



1 ig^_ f d 4 kYS(p - kYf V F p ^ 
2 CAjL (27r) 4 V " fc 2 (fc + g) 2 



regular part = (-c a )[——t\ / — — ^ (2.11) 



with ca the Casimir operator for the adjoint representation [ ca = N in SU(N) ] and one 
factor of 2 for the two pinching terms from Eq. (2.8) and Eq. (2.9). When we add to the 
usual propagator graphs (in the Feynman gauge) the pinch contributions of Eq. (2.10), 
together with an equal contribution coming from the mirror-image graph of l.b, we find 
the new gauge invariant self-energy Yl^ u (q), given by 

*W?) = (q 2 9^ - q^Mq) (2.12) 

where 



Tl( q ) = -bg 2 M^ r ) (2.13) 



and b = ^3 is the coefficient in front of — g 3 in the usual one loop /3-function. We see 
that the modified propagator has a gauge independent self-energy and only a trivial gauge 
dependence originating from the tree part given by Eq. (2.1). After all pinch contributions 
from graph l.b have been allotted to the new gluon self-energy, the rest of graph l.b, 
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namely the expression in Eq. (2.11), is genuinely vertex-like and must be added to the 
usual QED-like graphs (Fig. 2a). The final gauge-invariant vertex is given by 

f« = *MV^ [ d 4 kys( P -k)rr^ CA Nd f r d 4 kys( P + q - k)^s( P - k) ll 

" (27r) 4U 2 J i k 2 (k + q) 2 1 2 c / J k 2 

(2.14) 

where r a is the fermion representation matrix, df its dimension, and Cf its Dynkin index. 
If we now act on the f 1 ^ given in Eq. (2.14) with q a and use Eq. (2.7) it immediately 
follows that r£ satisfies the following QED-like Ward identity: 

q^ = r a [E(p)-E(p)] (2.15) 

with E the usual quark self-energy in the Feynman gauge and p = p + q. [19] The 
application of the P.T. in a theory with spontaneous symmetry breaking like the Standard 
Model is significantly more involved [13,14]. The main reason is that the massive charged 
vector mesons couple to fermions with different masses, and therefore the elementary Ward 
identity of Eq. (2.3) gets modified to: 

ft = {jf> + ft - mi) — {i> — m 2 ) + (mi - m 2 ) 

(2.16) 

=S 1 (p + k)-S 1 (p) + (mi - m 2 ) 
The first two terms of Eq. (2.16) will pinch and vanish on shell, respectively, as they did 
before. But in addition a term proportional to mi — m 2 is left over. However, as it was 
shown in [13], any additional gauge-dependent contributions coming from Eq. (2.16) cancel 
exactly against contributions coming from graphs involving charged Goldstone bosons, 
whose couplings to fermions are also proportional to mi — mi- 

Finally, we conclude this section with a brief presentation of an alternative formula- 
tion of the PT introduced in [14] in the context of the Standard Model. In this approach 
the interaction of gauge bosons with external fermions is described in terms of current 
correlation functions, i.e. matrix elements of Fourier transforms of time-ordered products 
of current operators [20]. This is particularly economical because these amplitudes au- 
tomatically include several closely related Feynman diagrams. When one of the current 
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operators is contracted with the appropriate four-momentum, a Ward identity is triggered. 
The pinch part is then identified with the contributions involving the equal-time commu- 
tators in the Ward identities, and therefore involve amplitudes in which the number of 
current operators has been decreased by one or more. A basic ingredient in this formula- 
tion are the following equal-time commutators, some of which we will also employ later in 
section 3: 

S(x - yo)[J^(x), J£(y)} = c 2 J^(x)5 4 (x - y) (2.17) 

S(x - yo)[J$r(x), Jti(y)} = -Jg(x)8*(x - y) (2.18) 

S(x - yo)[J$v(x), J£(y)] = J^(x)S 4 (x - y) (2.19) 
with J% = 2(J% + s 2 J^). On the other hand 

8(x -y )[4 i (x),j£.(y)} = (2.20) 

where Vi,Vj G {7, -Z} To demonstrate the method with an example, consider the vertex 
shown in Fig. lb, where now the gauge particles in the loop are W s instead of gluons 
and the incoming and outgoing fermions are massless. It can be written as follows (with 

r, = J ^r^( q ,k,-k-q) J d 4 xe ik *<f\T*[J#(x)J^(0))\i> (2-21) 

When an appropriate momentum, say k a , from the vertex is pushed into the integral 
over dx, it gets transformed into a covariant derivative acting on the time ordered 
product < f\T*[jfy (x) «/^(0)]|z >. After using current conservation and differentiating 
with respect to the 6>-function terms, implicit in the definition of the T* product, we 
end up with the left-hand side of Eq. (2.18). So, the contribution of each such term is 
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proportional to the matrix element of a single current operator, namely < f\ J%\i >; this is 

precisely the pinch part. Calling the total pinch contribution from the of Eq. (2.21), 
we find that 

rj = -g s d ww (Q 2 ) < f\4\i > (2.22) 

where 

/rl 4 k 1 

(tA) o 2 T" ( 2 -23) 

Obviously, the integral in Eq. (2.23) is the generalization of the QCD expression Eq. (2.8) 
to the case of massive gauge bosons. 



3. Gauge dependence of the conventional MDM form factor. 

In this section we explicitly show that the MDM form factor extracted from the 
conventional vertex graphs, calculated in a general renormalizable (R^) gauge, depends on 
the gauge-fixing parameter £, and is therefore not suitable for comparison with experiment. 
We then isolate the gauge dependent contribution to the MDM originating from the vertex- 
like pinch parts of box graphs, calculated in the same class of gauges. It turns out that 
the latter gauge-dependent contributions to the MDM cancel exactly against the gauge- 
dependent parts of the conventional vertex graphs, thus giving rise to a gauge independent 
expression for the top MDM. It is important to emphasize that even though it has been 
proved [14] that the inclusion of the vertex-like pinch parts is necessary in order to define 
a gauge independent jtt vertex, it has not been shown explicitly that the conventionally 
defined expression for the MDM is indeed gauge dependent. The obvious alternative (which 
we will show not to be true) is that, even though the usually defined vertex has an overall 
gauge dependence, its part proportional to cr^ v q u could be fortuitously gauge independent. 
In addition, the explicit computation of the gauge dependence of the conventionally defined 
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MDM enables one to appreciate quantitatively the significance of the issue at hand. The 
tree-level vector-meson propagator A flu (k) t in a general covariant gauge is given by 

Ai.„ W = [^- (3-D 

with i = w,z, 7, and M 7 = 0. The propagator D c (k) of the two charged Goldstone bosons 
is given by 

°^ = W^Ml ^ 
and the propagator D n (k) of the neutral Goldstone boson is 

b " w = s^W (a3) 

The propagator Djj(k) of the physical Higgs particle is gauge-independent at tree-level 



and therefore the graph of Fig. 3. 16 has no gauge-dependent contribution. Finally, we 
denote the quark propagators by 

S,(k) = y-L- (3.5) 

with j = t, 6, for top and bottom quarks. We further define c = cosOw and s = sinOw, 
where 9w is the Weinberg mixing angle, k = \Utb\ for the modulus of the C.K.M. mixing 
between top and bottom [21], and Cz = 1 — §s 2 — 75- Finally, the left-right projectors 
are given by P LjR = 1 =F 75 

We denote by the ^-dependent parts of the graphs shown in Fig. 3, with 
i = 1,2, ..,16 

We use the identity 

(e-i) 1 



["□> T7z\ ( 3 - 6 ) 



(k 2 - M 2 )(k 2 - £M 2 ) M 2V k 2 -iM 2 k 2 -M 2 ' 
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in order to isolate the ^-dependent parts of the graphs. Gauge-dependent terms propor- 
tional to are set to zero, when considered "sandwiched" between the on-shell target 
electrons. In what follows we omit the loop momentum integration symbol (^75) 2 / f^§ 5 
and we use dimensional regularization whenever necessary. We define the scalar quantities 
Di(£i), D 2 (Ci), and £> 3 (£i) as follows 

1 



2 



[(/ c + 9 )2_ M 2 ](A .2_ e . M 2 ) (3-7) 

D 3 m = 1 



[(k + qF-tiMfUV-tiM?) 

Clearly, £> 2 (£ = 1) = £> 3 (& = 1). 

The most involved graph is of Fig. 3, given by 

r« = KT^ TX A T P(k + q)/\ Xa {k) lp P L S b {k+p 2 ) la P L (3.8) 

Its ^-dependent part V p (£ w ) is given by 

Vf(t v ,) = <?\V} p + V2 P + V£ p ] (3.9) 

with 

t ,ip M , \ ^ p P L S b {k +P2)[m t PR - m b P L ] 

V = «(i - 5* 



v so,, + g)2 _ ^ [{k + g)2 _ {k2 _ M2) 



V 2P _ „ n t \ i m b p R ~ m t P L ]S b (k + p 2 h^P L , . 

" [ Kw} [{k + q y - miw - Mi){k* - t w Mi) ^ iUJ 

t/3P _ t x2 k^[m b P R - m t P L ]S b (k + p 2 )[m t P R - m b P L ] 

V — K{1 — l; w ) 



[(k + q y - Mi\ [(k + q y - i w Miw - Mi)(k* - $ w Mi) 

The superscript P in the formulas above stands for "pinch"; it hints to the fact that these 
contributions will eventually cancel against pinch parts coming from box graphs (see next 
section). It is important to notice the presence of the factor q 2 in the r.h.s. of Eq. (3.9). 
Clearly, Vf(£ w = 1) = 0. 
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The ^-dependent contributions vffi from each individual graph are: 



v fJ- 
v fJ- 



vf{t w )-vw-vp 

^P L S b (k+p 2 )[m t P R - m b P L ]D 2 ^ w ) - 
[m b P R - m t P L }S b (k + p 2 )^P L D 2 (£ w ) - \vj?> 



K 

~M, 



2k 



j^k^[m b P R - m t P L ]S b (k + p 2 )[m t P R - m b P L ]D 3 (£ w ) 



y(5) = J1_ 
» 3M.2 



ln p LSb{k + pi)[m t P R - m b P L ] - [m b P fl - m t PL\S b {k + p 2 )^^P L 



+ 2 llx P L D 1 {i w ) 



(3.11) 



^ (6) = df^ 7 ^ ~ 2m ^ C zSt(k +Pih 5 + 2m a5 S t (k + p 2 h II C z }D 1 (U - V™ 

y(7) = _ t/ (3) 

= ^(^)[m b P R -m t P L ]S t (k+p 2 )^S t (k+p 1 )[m t P R -m b P L ]D 1 ^ w ) 

V™ = l(^h5St(k+P2)l„S t (k + p 1 ) l5 D 1 (Z z ) 
V,^ = 



(i) 

The analogous ^-dependent contributions from the self-energy graphs T,\ of the external 
top quarks (Fig. (3)) are: 



K 



3M2 



ln p L S b{k + Pi)[m t P R - m b P L ] - [m b P R - m t PL]S b {k + p 2 )^^P L 



+ 2 7m P l D 1 (£ w ) 



_ V {14) 



V i 12) = ~^72^ C z - 2rn t ^C z S t (k + Pl ) l5 + 2m tTo S t (k + p 2 )^C z ]D 1 ^ z ) - 



y(13) = 



3M V 

4 3 1 
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(14) = 



ltxS t {pi)[m b P R - m t P L ]S t (k + Pi)[m t P R - m b P L ] + 

(3.12) 



+ [m b P R - m t PL}S t (k + p 2 )[m t PR - m b P L ]S t {p2)l^ 



£>i(U 



^(16) = o 

Adding the above equations together we see that Yllti Vi i] = Vf, so that 

^ ] =^\^=i) + V, P (U (3.13) 

In the equation above r^|(^ =1 ) represents the "regular" (hence the superscript R), purely 
vertex-like contributions in the Feynman gauge (£ w = 1). The term "regular" or "purely 
vertex-like" refers to what is left from the vertex graphs, after the propagator-like pieces 
have been removed through pinching, (see Eq. (2.11) and Eq. (2.14)) 

The calculation of the gauge dependence of the Zti vertex proceeds in an analogous 
way. The final answer is exactly the same as in Eq. (3.13), except that one must replace 
V£ by V^ z defined as 

V?* = (q 2 ~ M 2 )[V} P + V 2 / + V* p ], (3.14) 

e.g. 

^ w) =^\(^i) + V, Pz (U (3.15) 

It is important to notice the particular form of the gauge-dependent terms V/f and 
V^ z , namely that they are explicitly multiplied by q 2 and q 2 — M 2 , respectively. As we will 
see in the next section, this is precisely the characteristic structure of vertex-like contri- 
butions originating from box graphs, after pinching. Moreover, the gauge-dependent parts 
vanish, when the incoming gauge bosons are on-shell (q 2 = and q 2 — M 2 , respectively), 
as they should. 
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We next proceed to extract from Eq. (3.13) the form factor F(q 2 ) proportional to 
o'nuQ 1 ' '• We use the Gordon decomposition 

to convert terms proportional to {p\ + p 2 )^ into linear combinations proportional to 7 M 
and a^ v q u . The final answer is: 

F 2 (q 2 )=F 2 (q 2 )\ { ^ =l) + f 2 (Uq 2 ) (3.17) 

where -^2 (q 2 ) I (^^=1) originates from the first term in the r.h.s. of Eq. (3.13), whereas 
f2(iw,(l 2 ) originates from the second one, and has the following closed form: 



f(c 2 2q 2 m 2 f\ r 1 -* 

j2{€w,q ) = M2 dx dy 



xR(l, U + yR(U l)-(x + y)R(l, 1)+ 



(3. 
with 

R(Zi, 6) = [(x + yfm 2 - q 2 xy - (x + y)(m 2 + m\) + x^M 2 + y&M 2 + m 2 b ]~ l (3.19) 

Clearly, F 2 (q 2 ) of Eq. (3.17) is ^-dependent, through the explicit dependence of f 2 (i w , <? 2 ) 
on £ w . This concludes the proof that the conventionally defined MDM form factor is 
gauge-dependent, for both an incoming photon or Z, and therefore unsuitable in general 
for comparison with experiment. 

It is interesting to mention that there is no contribution from either parts of Eq. (3.13) 
proportional to the CP violating term 07^(^75. Therefore, the usual statement that there 
is no one-loop contribution to the EDM is indeed a gauge-independent one. 



8) 
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4. Gauge-independent MDM form factors via the PT 



In this section we show how the use of the PT can give rise to gauge-independent 
expressions for the MDM form factors, for both an incoming photon and Z. The important 
point is to recognize that the box-like graphs, like the one shown in Fig. 2b, contain gauge- 
dependent contributions which are kinematically equivalent to a 'jti or a Ztt vertex. When 
all these contributions are added to the conventional vertex graphs, the gauge dependence 
cancels exactly. 

The only kinematically relevant box graph is the one with an intermediate W + W~ 
pair, shown in Fig. 2b . [22]. Its ^-dependent contributions B are given by [23] 

B(U = 9 2 e-^P L e + g^[V^ + V 2P + V* p ] + ... (4.1) 

The dots in Eq. (4.1) indicate terms that belong to one of the following categories: 

a. Propagator-like, to be allotted to the 77, 7Z Z7, and ZZ self-energies. 

b. Vertex-like of the form g 2 [V* p + V 2P + V* p ]g v H lljL P L t; the V,} p , V 2P and V* p are 
given in Eq. (3.10) after replacing t <-> e and b «-> v. These contribution can be combined 
with the conventional , ye + e~ and Ze + e~ vertex graphs to construct gauge-independent 
■ye + e~ and Ze + e~ vertices. 

c. Purely box-like, to be combined with the rest of the box graphs. 

Clearly, none of the categories listed above contributes kinematically to the MDM 
form factors. 

We now write the expression given in Eq. (4.1) as a linear combination of a r ye + e~ 
and a Ze + e~ bare vertex, using the simple identity 

2 2 
9 -e + ltl P L e- = g 2 s 2 e+ lfl e- + ^e+^[P L - 4s 2 }e~ (4.2) 
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so that B becomes (we omit the external e + and e spinors) 

B = - ^V« P + V 2 / + V* p )] - { ^]^ L M2 4S2] [if ~ Ml){Vl P + v£ p + V* p )] 

q q ivi-z 

g 2 s 2 r v p _ d)r[PL-^} 

q 2 " q 2 -M 2 " 

(4.3) 

If we now add the first term on the r.h.s. of Eq. (4.3) to Eq. (3.13), and the second term 
to Eq. (3.15), we find the gauge-independent ^yti and Zti vertices, and A M , respectively. 
They are simply 

rM = rJ?| (c .=i) (4.4) 

and 

A, = A?I(C.=D (4-5) 
thus, the gauge-independent MDM form factor F(q 2 ) is given by: 

F(q 2 ) = F(q 2 )\, w=1 (4.6) 

for both incoming photon and Z. 

The gauge-independent •yti and Zti vertices given in Eq. (4.4) and Eq. (4.5) coincide 
with the the conventional vertex graphs of Fig. 3, calculated in the £ w = 1 gauge, but with 
the propagator-like pieces of graph 3.1 removed through pinching. In particular, the gauge- 
independent MDM form factors can be obtained if one just calculates the conventional 
MDM form factors in that same gauge. This results may come as no surprise to the 
reader familiar with the PT. Indeed, as we explained in section 2, in calculating gauge- 
independent Green's functions (a vertex in our case) any gauge choice is as good as any 
other, as long as one identifies and judiciously allots the pinch contributions. If one decides 
to calculate the gauge-independent vertex in the £ w = 1 gauge, the only possibility for 
pinching comes from the three-boson vertex of graph 3.1, which gives rise to a propagator- 
like piece proportional to 7 m -Pl- In particular, there are no pinch contributions (propagator 

17 



or vertex-like) coming from box graphs, since there are no longitudinal terms in the bare 
gauge boson propagators appearing in the boxes. Since the only term removed from the 
vertex has no part proportional to c^q", and hence no contribution to the MDM, the 
answer for the gauge-independent MDM is obviously the one recorded in Eq. (4.6). The 
fact that the gauge-independent answer in the MDM case coincides with the result in the 
Feynman gauge is in a sense fortuitous, and is certainly not true in general. Indeed, several 
other gauge-independent modified Green's functions, constructed via the PT, turned out 
to be different from the corresponding result in the Feynman gauge, or any other popular 
gauge for that matter. 

Finally, we conclude this section with a comparison of the gauge-independent = 1) 
answer for the MDM, to the answers computed for different values of £ w , as given in 
Eq. (3.17) and Eq. (3.18). We chose m t = 130 GeV and M H = 150 GeV. The results 
shown in Fig. 4 are rather impressive and convincingly demonstrate the importance of the 
issue at hand. Indeed, the residual gauge dependence is not just a theoretical nuisance, 
but has a huge impact on the predictions of the theory, affecting them both quantitatively 
and qualitatively. In Fig. 3 the real part of the MDM form factor for an incoming photon 
is plotted as a function of y/s. The curve marked (a) is the gauge-independent answer, 
(b) corresponds to the gauge choice £, w = 10, (c) to £ w = 100, (d) to t; w = 1000, and 
(e) to the unitary gauge (£ w — > oo) [24]. We see how the curves start approaching the 
unitary limit (e) , as £ w becomes large. The gauge-dependent answers are typically at least 
one order of magnitude larger than the gauge-independent answer of (a). The reason for 

2 

this artificial enhancement is the presence of the factor -Xr in front of the expressions in 

w 

Eq. (3.18), which is huge due to the massiveness of the top quark. Indeed, already at the 

2 2 

threshold of top pair production, for m t = 130 GeV, we have that tW- = 4-^ ps 10.5. If 

w w 

instead of tt we considered bb production the corresponding factor would be of the order of 
10 -2 . The phenomenological implications of such an enhancement are clear. For example, 
any dependence of the final answer on the mass of the Higgs boson gets totally altered. 
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Because of the same reason, the asymptotic behavior of the MDM form factor in a generic 
gauge is totally different than that of the gauge independent case; as we see from Fig. 4, 
the MDM curves increase monotonically as y/s — > oo, whereas, on the contrary, the gauge- 
independent curve approaches zero. It is interesting to notice that within the entire class of 
gauges, only the gauge-independent answer decreases asymptotically [25]. In addition, 
as we can see in curves (b) and (c), extra ^-dependent thresholds are introduced, which are 
artifacts and do not correspond to anything physical. Even though it is very unlikely that 
anyone would calculate in any of the finite gauges (£ = 10, 100, 1000) considered above, it 
is important to realize that the popular unitary gauge (e) yields a totally wrong answer. 
The same is true for another popular gauge, the Landau gauge, with £ w = (not shown). 
There, in addition to the usual artifacts, one has extra infrared divergences, due to the 
presence of massless poles. 



5. The Ward identity for the gauge-independent ytt vertex 

In this section we prove a QED-like Ward identity, which relates the gauge-independent| 

7ti vertex to the top quark self-energy E t in the Feynman gauge. The Ward identity is 

? Ai f At ( Pl ,p 2 ) = ^[E t (p 2 )-E t (p 1 )] (5.1) 

which is the one-loop generalization of the trivial tree-level result QtQ^^n = Qt[(jfe —nT>t) — 
(jfi — wt)], with Q t = g the electric charge of the top quark. It is important to emphasize 
that it is only after the propagator-like contributions of graph (3.1) have been removed 
that the validity of Eq. (5.1) becomes possible. To see that, let us recall the decomposition 
of the tree-level r yW + W~ vertex, given by Eq. (2.5) and Eq. (2.6). The graph , in the 
Feynman gauge, is (we again omit the (7^75) 2 / f^l factor in front): 

r?> = + O^^T^[ (fc2 _ M i )[{ l +q) > _ Mi] ] (5-2) 
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The part of Eq. (5.2) proportional to , after pinching out the internal quark propa- 
gator, gives rise to propagator-like contributions, which, as explained in section 2, will be 
allotted to the self-energies of the incoming neutral gauge bosons (77, ZZ, 7Z, Zj). How- 
ever, unlike the QCD residual contribution survives, because the internal (6) 

and external (£) quarks have different masses. This extra contribution is 

v r _ {[m t P R - m b P L ]S b {k)inP L + 'j fl P L S b (k)[m b P R - m t P L ]} 

^-Ml)[(k + q y-Ml] (5.3) 

= -(r^ + r( 3 )) 

and cancels entirely against the contributions of graphs Ip and Y^' . The above obser- 
vation is crucial for deriving Eq. (5.1). So, the purely vertex-like contribution Y^ L ' from 
is given by 

t (vl) _ K ^ P rP L s b (k)yp L _ 3 

After these considerations, proving the Ward identity is rather straightforward. In 
that vein, it is far more convenient to act with on each graph individually, instead of 
first calculating them, and then acting with on the final answer. Using Eq. (2.7) we have: 

<[r (^) +r ( 2)+r (3) ]=E (l) (p2) _ E (l) (pi) 

^r( 5 ) = -i[Ej 1) ( P2 )-E; i) ( Pl )] 

q^ = 2 -[^\ P2 )-^\ Pl )] 

frp-l^M-^M] (5 ' 5) 



3' 

^r(f) = 2 [E (5) (p2) _ E (5 )(pi)] 
? Mrao) = 2p(6) (p2) _ E (6 )(pi)] 

Adding the left and right hand sides of Eq. (5.5) we arrive at the advertised result. Clearly, 
pinching is instrumental for the validity of Eq. (5.1). 
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6. Calculations and results 



In this section we report the results of the numerical evaluation of the integrals listed 
in the Appendix. These integrals correspond to the vertex graphs of Fig. 3, calculated in 
the £ = 1 gauge, which as we explained in section 4, coincides with the gauge-independent 
answer in the PT framework. In arriving at these expressions we used the on-shell con- 
ditions 1/1=1/2 = rn t (or equivallently pf = p| = m t)i an d the Gordon decomposition of 
Eq. (3.16). 

We have used very accurate integration routines contained in the Mathematica pack- 
age. Such routines can handle the singularities on the integration path which appear in the 
formulas. The stability of the numerical results has been tested against different choices 
of the parameters of the algorithm. 

In an attempt to stay close to the notation of [2], we denote by Re(C 7 ) and Im(C 7 ) 
the real and imaginary parts of the MDM form factor for an incoming photon. Similarly, 
Re(C z ) and lm(C z ) denote the respective quantities for an incoming Z. The results are 
shown in Fig. 5, Fig. 6, and Fig. 7. In what follows we always assume that the aforementioned 
quantities can be individually extracted from e + e~ experiments, through the methods 
discussed in detail in [2]. In all our calculations m t > 130 GeV ', so the top is expected 
to decay before it can form a bound state [26], and therefore non-perturbative QCD 
effects should be numerically not very important. Even so, the one-loop QCD vertex 
corrections (graph 7 in Fig. 3 with 7 <-> gluon) dominate in general the contribution of all 
other graphs. This dominance is due to the fact that the strong coupling a s , although it 
is small enough for perturbation theory to be trusted, it is still significantly larger than 
both the electro-magnetic and the weak couplings, which multiply all other graphs. We 
used a s = 0.1, a = cos9 w = 0.885, M w = 80.6 GeV, M z = 91.1 GeV and 

nib = 4.5 GeV. In addition, in the computations we report, we always place ourselves 
sufficiently above the respective ti threshold, so that non-perturbative threshold effects 
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may also be omitted [27]. Of course, since the mass of the top is a free parameter which 
we vary in our computations, the location of the threshold (s t hr — 4m|) is also variable. 

In Fig.5, Re(C 7 ), Im(C 7 ), Re(C z ), and Im(C z ) are plotted as function of y/s for 
m t = 130 GeV and two different values for the mass of the Higgs: a light one (run = 
150GeV) and a heavy one (m H = bOOGeV). In Fig.6, Re(C 7 ), Im(C 7 ), Re(C z ), and 
Im(C z ) are plotted as function of -Js with a heavier choice for the top mass, namely mt = 
200 GeV, and the same values for the Higgs masses as in Fig.5. Finally, in Fig. 7 Re(C 7 ), 
Im(C 7 ), Re(C z ), and Im(C z ) are plotted as a function of m#, for two different values of 
the top mass, m t = 130 GeV and m t = 200 GeV. We fixed the value of y/s at 0.5 TeV , 
exactly as in the analysis of [2]. It is important to notice that in all cases shown in Fig.5 
and Fig.6 the form factors approach zero, as increases, at least for the range of the 
parameters that we have explored. 

As we see in Fig.5, both Re(C 7 ) and Im(C 7 ) are smooth functions of slightly 
above y/s = 300 GeV for both values of run- On the other hand both Re(C z ) and Im(C z ) 
display a spike (peak) at around yfs = 600 GeV when mjj = 500 GeV, whereas they are 
smooth and have no spikes when m# = 150 GeV. Exactly the same qualitative behavior 
is also displayed by the graphs of Fig.6; the only difference is an overall enhancement in 
the values of all quantities shown, since the top is now heavier. Finally, in Fig. 7 both 
Re(C 7 ) and Im(C 7 ) are smooth functions of m#; on the contrary, Re(C z ) and Im(C z ) 
are rather sensitive to the values of ran, especially if the mass of the top is large. 

The Feynman graph responsible for the behavior described above is the non-Abelian 
graph denoted by "special" in Fig. 3 (and its mirror graph, not shown), which exists only 
if the incoming gauge boson is a Z. On the other hand, the Abelian graph in Fig 3.10, 
common to both incoming 7 and Z, is of no particular importance. The importance of the 
"special" graph is due to the appearance of a threshold effect, as soon as the momentum 
y/s, carried by the incoming Z, reaches the critical value yfsl, given by yfs^ = M z + run, 
namely the threshold for Z + H production (e + e~ — > Z — > Z + H) [28]. It is now easy 
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to understand the behavior of Re(C z ) and Im(C z ) in Fig.5(c,d) and Fig.6(c,d): in the 
first case (m t = 130 GeV), for m# = 150 GeV we have yfs^ = 241 GeV ', which is below 
the ti production threshold y/s t hr = 260 GeF. So, the threshold effect of the "special" 
graph does not show up in the plots. On the contrary, for mu = 500 GeV, we have 
that yfsZ = 591 GeV, which gives the effect we observe in our plots, for both Re(C z ) and 
Im(C z ). Exactly the same happens in the second case (m t = 200 GeV), since the value 
of yfsl = 591 GeV is still larger than y/s t hr = 400 GeV. Of course, the value of yfsZ does 
not depend on the choices of mt; what depends on the value of mt is whether or not the 
ifZ-threshold effect at y/s^ appears in the plots. 

The threshold effect seen in both Fig. 7c and Fig.7d is also due to the "special" graph. 
For a given value of y/s we expect a threshold effect as soon as mu = mnlthr = M z , 
thus, for incoming y/s = 500 GeV one expects an effect at mn\thr = 500 GeV — M z xs 
409 GeV; this is indeed what we observe in Fig. 7c and Fig.7d. In particular, the position 
of m 1 ^ remains unchanged as it should, for both values of mt we examined. We emphasize 
that due to such threshold effect, some of our gauge-invariant form factors are much more 
sensitive to the value of the Higgs mass than the cross section as a whole [18]. For example, 
Re(C z ) of Fig.7c, for m t = 200 GeV (dark curve), changes from Re(C z ) = -1.4 x 10" 3 
at m H = 200 GeV to Re(C z ) = -3.9 x 10" 3 at m H = 400 GeV, a change of w 280%. 

Based on what we said above about the Zif-threshold effect of the "special" graph, 
the following possibilities may come up, if one measured Re(C z ) and Im(C z ) in the 
future (s max is the maximum energy of the collider): 

a) If 2m t — M z < mu < y/Smax, a threshold effect should appear in both Re(C z ) 
and Im(C z ) when plotted as a function of y/s, at a position y/s^ = mu + M z . 

b) If mu < 2m t — M z or m# > y/Sm ax , then no threshold effect should appear in 
either Re(C z ) or Im(C z , in the context of the SM. The second case is particularly interest- 
ing, especially in view of the current theoretical bias toward a relatively " light" Higgs. For 
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example if ran < 150 GeV, as favored by the Minimal Supersymmetric Standard Model, 
one should not see any HZ threshold effects, unless it turns out that m t < 120 GeV [29]. 

c) If a threshold effect appears in one of Re(C z ) or Im(C z ) but not in the other, or if 
more than one threshold effect appears in either or both of them, this should be interpreted 
as a signal of physics outside the SM. 

It is important to keep in mind that if one calculates the form factors in a gauge- 
dependent way, i.e. without pinching, in some cases one gets wrong predictions of threshold 
effects, related to unphysical, gauge-dependent thresholds. In these cases, the difference 
between the gauge-independent and the gauge-dependent prediction becomes, of course, 
dramatic. 



7. Conclusions 

In this paper we addressed some of the theoretical issues involved in the computation 
of off-shell form factors in the context of the Standard Model. Such form factors can 
be extracted, at least in principle, from cross-sections in future e + e~ experiments. In 
particular, we explicitly demonstrated that the conventionally defined MDM form factors 
of the top quark are gauge-dependent in the class of gauges. Due to the massiveness of 
the top quark, this gauge dependence turned out to be numerically very strong. We used 
the S-matrix pinch technique to eliminate the gauge dependence, and defined a gauge- 
independent MDM form factor, suitable for comparison with experiment. We proved a 
simple QED-like Ward identity, satisfied by the gauge-invariant 'jti vertex costructed via 
the PT. Finally, we computed the gauge-independent MDM form factors and presented 
the results for both real and imaginary (absorptive) parts. A very interesting feature of 
such form factors is their high sensitivity to the value of the Higgs mass, which enters the 
formulas as a parameter. For this reason, assuming that the top quark will be discovered 
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and detected before the Higgs particle, the form factor we have analyzed in the present 
paper contains interesting information about the Higgs boson. 

Of particular interest for the study of CP violation will be the experimental measure- 
ment of EDM form factors in the future. As we mentioned in the introduction, there are no 
EDM contributions within the SM up to three loops, but several models predict non-zero 
EDM already at one-loop. We are currently investigating issues of gauge invariance in the 
context of some of those models, with focus on the EDM form factors, and will present 
our results elsewhere. 
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9. Appendix 



In this Appendix we list all the contributions to the magnetic form factor F(q 2 ) of 
the top quark. Using the standard Feynman parametrization formula 



ABC~ J dX J{] 



l-x 



dy- 



(9.1) 



[Ax + By + C(l -x-y)f 

and subsequently the convenient change of variables t = x + y, tz = x — y, we can cast 
the MDM contributions of the graphs in Fig. 3 (except the one denoted "special") in the 
following form 



V 

(0 



dt t 



dzL {i) (t,z) } (9.2) 
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where L^(t, z) = L(t, z; m\, 777-2, Q 2 )- 

In what follows we will use the short-hand notation L(mi, 7772). In order to identify 
the contributions to the magnetic form factor F J (q 2 ) and F z (q 2 ), which arise from the •ytt 
and the Ztt vertices respectively, the coefficients in Eq. (9.2) carry a vertex label V = 7, Z. 
Obviously the integrand functions (z) do not have such a label, as they do not depend 
on which gauge boson enters the vertex. In the following list, the contributions are listed 
in the same order as in Fig. 3 and we use again the abbreviations s = sinOw, c = cosOw 
and k — \U t b\ 2 - In addition, we define C = (1 — |s 2 ). 



7 _ orv_ 7 7 7 _ 

JV (l+2+3) - 4s 2 /i (l+2+3) - U ^(1+2+3) ~ 1 °(l+2+3) ~ 1 

L (1+2+3) = L(M 2 ,m 2 ) (9.3) 



TV 7 

iV (4) 



3« 



8 M 2 s< 



A (4) 



2m? 



R 7 

^(4) 



—3m? — m? 



/^7 
°(4) 



m 2 . + m 2 



3k 



8 M 2 s 2 



2s c 



A (4) 



2mi 



-3m 2 — m 2 



°(4) 



m h + m 



L (4) =L(M 2 ,m 2 ) 



(9.4) 



N h 

3k 
8 s 3 c 



4 Z 

A (5) 



4 7 



= -2 + -s' 



-3 



/^7 
C (5) 



S f 5 ) = 3 



2s z 



°(5) - 



"I + ^ 



L (5) =L(m 2 ,M 2 ) 



(9.5) 



/V 7 - 



4 s 2 c 2 



4 7 — 9 



B 1 --3-— s 4 + -s 2 
^(6) - d 9 + 3 
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64 s 3 c 3 



A z -- 



8 C Bf &) = -C(C 2 + 11) 
L {6) =L{m 2 ,M 2 z ) 



Cf 6) =£(£ 2 + 3) 



(9.6) 



N 1 --- 
iV (T)" g 



N m = 



3 sc 



Af 7) = 



B 1 - -1 

^(7) - 1 



B {7) = 



°(7) 



C (7) - 1 



'(7) 



L(m 2 ,0) 



(9.7) 



_/yT = 1 

W 8 Mis 2 



3k 



64 Af2 s 3 c 



4 Z 



A 1 



R 7 

^(8) 



4(l--s 2 )(m 2 -m 2 ) 



L (8) =L(mlM 2 ) 



C (8) = m b+ ™ 2 t 



Cf 8) = -4m 2 (l-^ 2 ) + j|m 2 S 2 | 

(9.8) 



mi 



TV 7 = 
( 9 ) " 4M 2 ,s 2 



4 7) - n 



S (9) = 



°(9) - 1 



3 m 2 



32 M2 S 3 C 



A (9) 



B (9) 



L (9) =L(mlM 2 ) 



°(9) 



(9.9) 



mi 



( 10 ) 4 M^s 2 



iV ( 10 ) " 32 M 2 S 3 C 



3 m? £ 



z _ 3 m 2 
iV ( 10 ) " 32 M^ S 3 C 



^(10) ~ 



l (io) 



A z 

A {1Q) 







B U = 2 



^(10) 



tdZ 
^(10) 



L m =L{m 2 u M 2 H ) 
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^(10) - _1 



°(io) 



°(io) 



-1 



-1 



(9.10) 



For the "special" graph (last one in Fig. 3) and its mirror graph (not shown) the original 
Feynman parametrization given in Eq 9.1 was maintained, since the change of variables 
(x, y) — > (t, z), used for all other graphs, is not convenient. 

Finally, the one-loop QCD contribution is obtained from Eq. 9.6, by multiplying by 
the factor where a s is the QCD coupling at u. = M z , namely 

2a (9.11) 



QCD l 2a \ (7) 
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29. Of course in such a case the usual arguments about the non-hadronization of the 
produced tops will be invalidated, and non-perturbative QCD effects may become 
important. 

11. Figure Captions 

1) Graphs (a)-(c) are some of the contributions to the S-matrix T. Graphs (e) and 
(f) are pinch parts, which, when added to the usual self-energy graphs (d), give rise to a 
gauge-independent effective self-energy. 

2a) The two graphs defining the gauge-invariant vertex T^. The circle in the first 
graph indicates the part of the three-gluon vertex. 

2b) Vertex-like pinch contributions from a box graph, in a general gauge, with 

3) All the graphs contributing to the ^ti and Ztt vertices. 

4) The top quark MDM as a function of y/s, computed for several different values 
of the gauge-fixing parameter £ w (b,c,d,e), to be contrasted to the gauge-invariant answer 
shown in (a). The correspondance is as follows: £ TO = 10 (b), £ TO = 100 (c), £ w = 1000 (d), 
£ w -> oo (e). 

5) The case of a "light" top (m t = 130 GeV), for two different choices of the Higgs 
boson mass: m# = 150 GeV, and run = 500 GeV (darker curve). In particular: 

5a) Real part of MDM for an incoming 7, 

5b) Imaginary part of MDM for an incoming 7, 

5c) Real part of MDM for an incoming Z, 

5d) Imaginary part of MDM for an incoming Z. 



6) The case of a "heavy" top (mt 
boson mass as in Fig. 5 : 



= 200 GeV), for the same two choices of the Higgs 
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6a) Real part of MDM for an incoming 7, 

6b) Imaginary part of MDM for an incoming 7, 

6c) Real part of MDM for an incoming Z, 

6d) Imaginary part of MDM for an incoming Z. 

7) The dependence of MDM on the mass tjih of the Higgs boson, for a " light" and a 
heavy" top, at yfs = 500 GeV: 

7a) Real part of MDM for an incoming 7, 

7b) Imaginary part of MDM for an incoming 7, 

7c) Real part of MDM for an incoming Z, 

7d) Imaginary part of MDM for an incoming Z. 



32 



Edinburgh U. 93/530 
NYU-TH-93/05/02 
hep-ph/93 11284 
November 1993 



On 

> 



Gauge invariant top quark form factors 
(N ! from e + e experiments. 



> 

oo ■ J- Papavassiliou 1 and C. Parrinello 2 . 

1 Department of Physics, New York University, 4 Washington Place, 
New York, NY 10003, USA. 



Oh 



X 



2 Department of Physics, University of Edinburgh, Mayfield Road 
Edinburgh EH9 3JZ, Scotland, UK 



ABSTRACT 



Motivated by the possibility of experimental determination of top quark form 
factors in upcoming e + e~ experiments, as recently discussed by Atwood and 
Soni, we show at the one loop level that the conventionally defined magnetic 
dipole moment (MDM) form factor of the top quark is gauge dependent, in 
the class of renormalizable (R%) gauges. We explicitly calculate its gauge 
dependence, which, due to the massiveness of the top quark, turns out to be 
numerically very sizable. We show how to use the S-matrix pinch technique 
in order to define gauge independent form factors. The real and imaginary 
parts of the gauge independent MDM form factor are calculated and their 
dependence on the top quark and Higgs boson masses is discussed. Most 
noticeably, the dependence on the Higgs mass turns out to be stronger than 
that of the cross section as a whole. 
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